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Abstract 

We analyze the dynamics of a four-dimensional null hypersurface in a five-dimensional bulk 
spacetime with Einstein- Yang-Mills fields. In an appropriate ansatz, the projection of the field 
equations onto the hypersurface takes the form of conservation laws for relativistic hydrodynamics 
with global non-abelian charges. A Chern-Simons term in the bulk action corresponds to anomalies 
in the global charges, with a vorticity term arising in the hydrodynamics. We derive the entropy 
current and obtain unique expressions for some of the leading-order transport coefficients (in the 
abelian case, all of them) for arbitrary equations of state. As a special case and a concrete example, 
we discuss the event horizon of a boosted Einstein- Yang-Mills black brane in an asymptotically 
Anti-de-Sitter spacetime. The evolution equations in that case describe the hydrodynamic limit of 
a conformal field theory with anomalous global non-abelian charges on the Anti-de-Sitter boundary. 

PACS numbers: 04.70.-s, 47.10. ad, 11.25.Tq 
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I. INTRODUCTION AND SUMMARY 



The hydrodynamics of relativistic conformal field theories (CFTs) has attracted much 
attention, largely in view of the AdS/CFT correspondence between gravitational theories 
on asymptotically Anti-de-Sitter (AdS) spaces and CFTs jlj] (for a review see Q). Hydro- 
dynamics is an effective description of the long distance field theory dynamics and applies 
under the condition that the correlation length of the fluid l cor is much smaller than the 
characteristic scale L of variations of the macroscopic fields. The AdS / CFT correspondence 
implies that the long wavelength dynamics of gravity is dual to the CFT hydrodynamics. 
Indeed, it has been shown in [3[ that the (d+ l)-dimensional CFT hydrodynamics equations 
are the same as the equations describing the evolution of large scale perturbations of the 
(d + 2)-dimensional black brane. The derivation in j^] is similar to the derivation of hydrody- 
namics from the Boltzmann equation, where the thermal equilibrium solution is the boosted 
black brane 4J. The limit of non-relativistic macroscopic motions in CFT hydrodynamics 
leads to the non-relativistic incompressible Navier-Stokes (NS) equations {5,6], and the dual 
gravitational description is found by taking the non-relativistic limit of the geometry dual 
to the relativistic CFT hydrodynamics 6]. 

In p| the membrane paradigm formalism 0-0 and an expansion in powers of the Knud- 



sen number l cor /L have been used to show that the dynamics of a membrane defined by the 
event horizon of a black brane in asymptotically AdS space-time is described by the incom- 
pressible Navier-Stokes equations of non-relativistic fluids. This has been generalized in [ll| 
to the relativistic CFT hydrodynamics case. The starting point in the membrane paradigm 
framework is an equilibrium (d + 2)-dimensional solution containing a timelike Killing vec- 
tor field and a stationary {d + l)-dimensional causal horizon. This solution is associated 
with a thermal state at uniform temperature. When a hydrodynamic limit exists, we can 
expand the general inhomogeneous black brane in the neighborhood of the causal horizon 
in powers of Z cor /L, assuming that there is no singularity at the horizon. For a black brane 
in asymptotically AdS, the set of Einstein equations projected into the horizon surface at 
lowest orders in l CO r/L is equivalent to the (d+l)-dimensional relativistic CFT Navier-Stokes 
equations. The membrane paradigm has previously been used as a tool for calculating trans- 
port coefficients 12Ml7| . The RG flow relation between the horizon coefficients and those of 



the boundary gauge theory have been recently discussed in 



18| 



It has been recently revealed that the hydrodynamics description exhibits an interesting 
effect when a global abelian symmetry current of the microscopic theory is anomalous. 



This has been first discovered in the context of the fluid/gravity correspondence |19l . |20 |. 
The Chern-Simons term in the gravity action, which corresponds to having an anomalous 
global abelian symmetry current in the dual gauge theory, has been shown to modify the 
hydrodynamic current by a term proportional to the vorticity of the fluid. At first sight the 
additional vorticity term seemed in contradiction with the second law of thermodynamics 



2l| . This, however, has been resolved by a redefinition of the entropy current in 22]. 



xperimental signals of the vorticity term from heavy ion collisions have been proposed in 
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24]. 



In this paper will analyze the dynamics of a four- dimensional null hypersurface in a five- 
dimensional bulk spacetime with Einstein- Yang-Mills fields. In an appropriate ansatz, the 
projection of the field equations onto the hypersurface takes the form of conservation laws 
for relativistic hydrodynamics with global non-abelian charges. A Chern-Simons term in 
the bulk action corresponds to anomalies in the global charges, with a vorticity term arising 
in the hydrodynamics. We will derive the entropy current and obtain unique expressions 
for some of the leading-order transport coefficients (in the abelian case, all of them) for 
arbitrary equations of state. As a special case and a concrete example, we will consider the 
event horizon of a boosted Einstein- Yang-Mills black brane in an asymptotically Anti-de- 
Sitter spacetime. The evolution equations in that case describe the hydrodynamic limit of 
a conformal field theory with anomalous global non-abelian charges on the Anti-de-Sitter 
boundary. 

We work perturbatively in the derivative /Knudsen number expansion, and write the 
equations to second order; equivalently, we consider the hydrodynamic currents and stress 
tensor to the first viscous order. We place all 4d fields on the null hypersurface, without 
explicit reference to the bulk spacetime or its boundary. When an AdS asymptotics is 
used for the unperturbed homogeneous brane solution, it provides the equation of state for 
the thermal system. Without solving the radial equations, we consider the most general 
inhomogeneous corrections to the horizon projections of the bulk fields. These corrections 
are then constrained by the requirement that the 4d constraint equations take the form 
of conservation laws. In the abelian case, the constitutive relations are obtained from this 
procedure in a unique closed form. For non-abelian charges, some of the terms are similarly 



obtained in closed form, with part of the conductivity matrix remaining ambiguous. 

The null surface formalism introduced here improves on the methods developed in [ll| 
by utilizing an intrinsic covariant divergence on the horizon. Our geometrical calculation 
is performed with the horizon's null normal as the basic vector variable. In viscous hydro- 
dynamics, this corresponds to the "frame" where the entropy velocity is the basic variable; 
however, the results are readily translated into the more standard Landau frame, which uses 
the energy velocity. 

The horizon dynamics equations define the viscous hydrodynamics of a relativistic CFT 
with non-abelian conserved currents and generalizes the known cases with conserved abelian 
currents. In the presence of a non-abelian Chern-Simons term, the horizon dynamics gives 
the hydrodynamics equations in the presence of anomalous non-abelian global symmetries. 
The result generalizes the hydrodynamics equations of field theories with triangle anomalies 
of abelian currents, discussed above. 

The thermodynamics of the AdS black brane is described by the equation of state 



T, chemical potentials \i a and pressure p are derived from the equation of state through the 
usual thermodynamic identities. 

In the present calculation, we will make no use of the particular equation of state ( 11.11) . 
Instead, we study a generalized ansatz for the horizon fields, corresponding to an arbitrary 
equation of state. Thus, in addition to the concrete case of the AdS black brane, we study 
more generally the structure of hydrodynamic systems which can be encoded holographically 
on a null hypersurface with Einstein- Yang-Mills fields and an Einstein- Yang-Mills- Chern- 
Simons action. In particular, do not use the conformal symmetry of the AdS asymptotics 
and of eq. (11.11) . 

We now proceed to a summary of the results, expressed in the Landau frame. For the 
viscous stress-energy density, we find the familiar expression 



where is the energy 4- velocity, and h^ v is the 4d background metric for the hydrodynamics. 
Our formalism is such that a curved h^ v can be allowed for with no additional complications. 
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(I.l) 



where e, s and n a are the energy, entropy and charge densities, respectively. The temperature 




(1.2) 



ix pu is the shear tensor, defined as 

V, = P^D {p u a) - l -P pv D p u p , (1.3) 

where P pv = h pu + u^u,, is the projector orthogonal to u M , and D p is the covariant derivative 
with respect to h pv . Eq. (II. 2[) implies the celebrated value 77 = for the shear viscosity, 
and a vanishing bulk viscosity £ = 0. 

For the (non-abelian) charge currents, we obtain the expression 



7 m _ ,/Z77 I „ „a« T l sl/3 ( Zn M + M^V^pFn^ 
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(1.4) 



where cj m = |e M ^ po T u v d p u a is the fluid's vorticity, /3 a fe c are the Chern-Simons/anomaly coeffi- 
cients, and Z ab is the projector onto the charge subalgebra Z(n) which commutes with the 
local charge density n a . The distinctively non-abelian contribution is concentrated in the 
coefficient matrix a^ b , whose indices lie in the orthogonal complement Z- L {ri) of Z(n). Its 
symmetric part cx^ is positive semi-definite, as required by the second law of thermody- 
namics. Other than that, our method places no restriction on the function a^ b (s,n c ). The 



constitutive relation (II. 4j) reduces to the known results [19|, [20| in the abelian case, where we 
substitute Z ab = 5 a t> and a^b = 0- The abelian case is also reproduced in the limit of small 
charges or weak couplings, as we will see in section I VI I 

We will find that the undetermined part a^ b of the conductivity matrix depends on the 
inhomogeneous corrections to the bulk gauge potential. In the AdS/CFT context, these are 
fixed by the radial part of the field equations together with the condition: 

lim A a — . (1.5) 

r—¥oo 

On the event horizon of the AdS black brane, the resulting corrections to A a ^ are given 
by some functional of the hydrodynamic fields. Without solving the radial equations, we 
cannot know the details of this functional. Outside the AdS context, our result is simply 
that the transport coefficient a^b can be arbitrary in Einstein- Yang-Mills-Chern-Simons hy- 
drodynamic systems (and at leading order, it is the only arbitrary transport coefficient). In 
the AdS / CFT case, we obtain algebraic relations between the corrections to the non-abelian 
gauge potential on the horizon and the conserved current J£. This current is given on 



the other hand by the AdS/CFT prescription as: 

J a = -— lim r 2 A a u . (1.6) 

2,71 r—^oo ^ 

Thus, we obtain algebraic relations between the corrections to the gauge potential on the 
horizon and on the boundary. We verified these relations numerically against the bulk 
first-order calculation for the £77(2) case carried out in 26]. 
Finally, our result for the entropy current reads: 

s V nysu + 2 4/ 37r(e + p) ^ e + p J U T) + 3vr(e + p) ' 1 j 
This result generalizes the entropy current derived previously on the gravity side of 
AdS/CFT using properties of the uncharged black brane solution [25]. 

The focusing equation of the null horizon requires the divergence of this entropy current 
to be non-negative, which translates into an entropy production rate: 

V = \d„S» = V^h (To ah P^Dj^Dj^ + ^v^) • (1-8) 

The Chern-Simons/anomaly terms in (11.41) and (11.71) reproduce the results of the general 
thermodynamic argument in 22], which was in turn motivated by the specific results from the 
AdS black brane. As mentioned above, our calculation yields these results for an arbitrary 
equation of state. Furthermore, we generalize the work in {22]] by considering hydrodynamics 
with non-abelian as well as abelian charges. 

The paper is organized as follows. In section [XT] we present the geometrical framework, 
we define the fields and field equations in the 5d bulk and their projections onto the event 
horizon. Section ITO1 describes the homogeneous (thermodynamic) black brane solution from 
which we begin and its various parameters, and outlines our approach to the inhomogeneous 
brane. In section HVl we derive the leading-order (ideal) hydrodynamic equations. The prin- 
cipal weight of the paper lies in section [V] where the second-order (viscous) hydrodynamic 
equations are derived. In section IVIl we discuss the results. In appendix |A] we provide a 
summary of the notation used in the paper. 



II. THE GEOMETRICAL FRAMEWORK 



We consider the Einstein- Yang-Mills fields: the metric g^B and the gauge potential A a A for 
a gauge group Q. These fields live in a 5d spacetime with a null horizon %. We consider the 
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relevant restrictions of the fields into the 4d geometry of the horizon, and develop the general 
formalism which will be applied to the hydrodynamic ansatz in the following sections. 

A. Horizon geometry 

In the bulk spacetime, we choose coordinates of the form x A = (r, The x M coordi- 
natize the horizon Ti; r is a transverse coordinate, with r = on %. d^r is a null covector 
tangent to the Ti. When raised with the metric, it gives a vector field £ A = g AB dsr which is 
both normal and tangent to "H, and tangent to its null generators. In components, we have 
£ A = (0,^ M ). These choices fix the following components of the inverse bulk metric on %: 

g rr = q. g r„ = t _ (jj -q 

The pullback of qab into "H is the degenerate horizon metric 7^. Its null directions are 
the generating light-rays of "H, i.e. r y tiv £ u = 0. The Lie derivative of 7^ along £ M gives us 
the shear /expansion tensor, or "second fundamental form": 

o^ = \cn^, {u.2) 

which has the properties: 

o,u = e v ,- i% = o. (11.3) 

(H) 

We can write a decomposition of Q^ v into a shear tensor o\ w } and an expansion coefficient 6: 

V = °js? + \^ ■ (n.4) 

An equivalent condition which we found useful in practice is 

9 = \{G- l Y u 6, u , (II.5) 

where (G~ l ) flL ' is the inverse of any matrix G^ v of the form: 

G^u = A T|U , - b^K; b^ ^ . (H.6) 

We introduced the superfluous scalar field A for later convenience: it will turn out that a 
matrix of the form (1II.6|) coincides at leading order with the metric h^ u of the hydrodynamic 
dual. 
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The expansion coefficient 6 is related to the horizon area density current S^, which is the 
analog of the bulk volume density \J—g. This current is a vector density tangent to the null 
generators, with magnitude defined by: 

C \tl\V\ lillVI ]f±3,V3 *J U ^VQV\V1VZ ■ \ LL - { ! 

Here, e^ upa is the invariant 4d Levi-Civita density with components ±1, and e^pa is the 
inverse density with components ±1. We have the collinearity relation S p = vi* 1 , where v 
(v = 4s) is a scalar density. Note that unlike the magnitude of is fixed uniquely by the 
horizon metric. In our adapted bulk coordinates, the value of the 4d density v(x fl ) equals 
that of the 5d density \J—g evaluated at the corresponding horizon point (O,^). Finally, 
the divergence of the area current S 1 * is related to 9 by: 

v0 = dpS^ . {11.8) 

Since 7 MI , is degenerate, one cannot use it to define an intrinsic connection on the null 
horizon, as could be done for spacelike or timelike hypersurfaces. The bulk spacetime's 
connection does induce a notion of parallel transport in 7i, but only along its null generators. 
This structure is not fully captured by 7^; instead, it is encoded by the extrinsic curvature, 
or 'Weingarten map' 0^", which is the horizon restriction of S7a^ B '- 

0/ = v M r . (11.9) 

For a non-null hypersurface, the extrinsic curvature at a point is independent of the induced 
metric. For null hypersurfaces, this is not so. Indeed, lowering an index on M ^ with 7^ 
(and losing information in the process), we get the shear /expansion tensor 9^: 

e/ 7p , = V • ( IL1 o) 

This expresses the compatibility of the parallel transport defined by with the horizon 
metric 7^. Contracting with yields the surface gravity k, which measures the non- 
affinity of 

9/^ = nt . (11.11) 

It follows from fill. 1 OH - (III . 1 1 j) that given an arbitrary G^ v of the form (III.6p . 0^ can be 
written as: 

e/ = \e w {G- 1 y v + C pt- C pt = k . (11.12) 
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The covector c p encodes the 4 degrees of freedom in Q p v which are independent of 7^. 
In the hydrodynamics, these degrees of freedom roughly correspond to the velocity and 
temperature fields. 

The parallel transport along null generators defined by Q p u can possess a Ricci-type 
curvature, described by a 4d covector density. In bulk terms, this quantity is given by the 
projection R pu S u of the bulk Ricci tensor Rab] therefore, it is coupled to matter by the 
Einstein equation. On the other hand, it is related to Q p u through a null version of the 
Gauss-Codazzi equations. We will now describe this relation, and bring it to a convenient 
form. Our discussion follows and expands on the formalism of 

In their usual form, the Gauss-Codazzi equations involve the divergence of the hypersur- 
face's extrinsic curvature. Our null horizon, however, does not possess an intrinsic connec- 
tion, and a covariant divergence of Q p u cannot be defined. The solution is to define a tensor 
density constructed out of Q p u in the following manner: 

Q/ = <©/-/^), (11.13) 

with the properties: 

l»Qv = 0; Q/7p, = Q/7p,u • (H.14) 

It turns out that the degenerate metric 7^ is sufficient to define the divergence of any 
density which satisfies AIL 1411 . We begin by defining a raised-index version of Q^, which 
satisfies: 

Q^ = Q^; i, P Q pu = Q, u - (11.15) 

The properties ( III. 151) do not define Q pu uniquely, but only up to multiples of l p l v ] this will 
be good enough for our purpose. The horizon-intrinsic covariant divergence of Q p can now 
be defined as follows: 

V„Q M " = d v Q» - \Q up d^ lvp . (11.16) 

Since the contraction of d fl 'y up with l v l p is zercQ, the potential ambiguity in its contraction 
with Q up drops out, and we have a well-defined expression. Using this divergence, the null 
Gauss-Codazzi equation can be written as: 

R^S" = V„Q M " - vd„9 . (11.17) 

1 This follows directly from the identity d^ up i v t p ) = 
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We found it convenient to use the decomposition (1II.12j) of Qp u when evaluating eq. 
(III. 171) . The corresponding decomposition of Qp reads: 



Q/ = XvO^G^r + M" - c p S%), (11.18) 

where each of the two terms separately satisfies eqs. (jII.14p . We can therefore consider their 
covariant divergences separately. Let's start with the first term. A raised-index version of 
\v6 fJiP (G~ 1 ) pu satisfying (1II.15j) can be obtained by simply raising with A(G~ 1 ) M ^; further- 
more, eq. (1H.3|) and the relation [G~ x Y v b v ~ £ M allow us to replace dp^ vp in fill. 161) with 
d^\~ x G vp ). We get: 

% (Xve^G- 1 )^) = Df> (Xv9, p (G^Y") + vdd^lnVx, (11.19) 

(G) 

where Dp is the covariant derivative associated with the "metric" Gp V . The divergence 
of the second term in (111.18j) is a simpler matter; an object of this type has a covariant 
divergence which depends only on S p and c M , with no need for the full jp U : 

V u (cpS v - c p S%) = Cpd u S u + 2S v d [v Cp ] . (11.20) 

Summing up, we have: 

Rp V S v = D<® (XvOppiG- 1 )^) + v9dp In + Cpd v S v + 2S v d [v c ll] - vd„6 . (11.21) 

The ^-component of Rp U S u can be expressed more directly, as the null focusing equation: 

Rp U tt = -tdpO + k9- -9 2 - \ 2 {G^ 1 Y P {G~ 1 Y a a { ppa i p^ ) . (11.22) 

3 

Once again, the shear-squared term could be defined more directly using the tangent bundle 
modulo displacements along light-rays. The form given here is more convenient for practical 
calculations. 



B. Yang-Mills 

The bulk gauge potential A\ is associated with a field strength F% B = 28[aA^ + 
f a b c A b A A c B . The structure constants f a b c = f[ a b c ] include the coupling strengths of each 
simple piece of the gauge group. The gauge group in the 5d spacetime induces a 4d gauge 
group on %. The corresponding gauge potential Ap[ is the pullback of the bulk potential 
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A\ into H. Likewise, the field strength F* v = 1d^A a u] + j a hc A\A c v associated with A a ^ is the 
pullback of F\ B . 

Another restriction of F AB into % with a well-defined 4d meaning is given by S v F av A . 
This quantity is tangent to the horizon. The upper index can therefore be restricted into 
the tangent bundle of "H, giving a 4d vector density. In our adapted coordinates, it equals 
y/—gF^. Inserting a numerical factor, we define: 



ja„ _ V_9 F ar^ = V t ( F a m + f « ( IL 23) 

An Air p 

where we used ( III. Ill in the second equality. At the ideal-fluid approximation, this will 
coincide with the hydrodynamic charge current. The numerical factor is chosen to make the 
chemical potential of the fluid agree with the "electric potential" t^A® on %. Note that j% 
is covariant under Ad gauge transformations on T-L. The two quantities A® and on the 
horizon play a role similar to that of the non-normalizable and normalizable parts of A a ^ on 
the AdS boundary. The quantities F^ v and j£ are not independent, and do not span all the 
components of F% B on H. This is a consequence of the null nature of %. Explicitly, we have 
the following relation, analogous to ( III. 1011 : 

l^r = ^F^ . (11.24) 

The single degree of freedom in j£ which is independent of corresponds to the hydrody- 
namic charge density. 



C. Field equations 

1. Bulk equations 

The bulk fields satisfy the Einstein equation with a cosmological constant and a Yang- 
Mills stress tensor: 

Rab ~ Aqab = StvTab , (H.25) 

Tab = ^{F a A cFaB C - \F a CD F™ g AB ), (11.26) 

coupled to the Yang-Mills equation with a Chern-Simons term: 

9b (V=9F a AB ) + V^9f ab c A B F c AB - P abc e ABCDE F BC F c DE = , (11.27) 
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where ^ is the invariant 5d Levi-Civita density with components ±1. The Chern- 

Simons coefficients f3 a b c = P( a bc) can be any gauge-invariant constant tensor. In particular, 
they can mix different U(l) and simple pieces of the gauge group Q. Examples include: 

• Any constant tensor for Abelian Q. 

• Pabc ~ for g = SU(N), N>3. 

• Pabc ~ ^(a^b^c) f° r S au S e g rou P G = G' x U(l) with G' simple, where is the index of 
the U(l) generator, and d' runs over the generators of G' ■ 

In general, the gauge-invariance of (3 a b c is encoded by the condition: 

f ed aPdbc + FWadc + r'cPabd = . (11.28) 



2. Horizon equations 



To restrict the Einstein equation (1II.25|) into H, we contract with S v and pull back the 



remaining index. The result is an equation between 4d covector densities. The ~ qab terms 
don't contribute, since qavS u ~ 9ab^ b = 9at, and that vanishes under pullback. We get: 

R^S" = 8n% u S u . (11.29) 

The LHS can be expressed in terms of lower-derivative geometric quantities on % via the 
Gauss-Codazzi equation (III. 211) . For the RHS, we have, using (IH.lj) and (III.23j) : 

v 

Air' 



V = 7z F ^au°r = f;ji . (n.30) 



Substituting into (III. 291) . we get the Einstein-Gauss-Codazzi equation: 

DP (Xve^G-y) + vOdp \nV\ + Cfl d u S u + 2S v d [vCtA - vd^O = 8irF«X . (11.31) 

The contraction of this equation with i^jv can be written more simply using the focusing 
equation ( III. 221) . We get: 

-1^0 + \e 2 - X^G^nG^ra^a^ = -£^ v f a . (11.32) 

O 1) 



The Yang-Mills- Chern- Simons equation ( II 1. 2 7 j) can be restricted to % by contraction with 



e A^po-e Mlypo "/24. Dividing further by 47r, we write the resulting equation as: 

+ ^cs) + f\cAUr + J%) + ^e^f b de A d AlF c pa = , (11.33) 
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where is defined in (III. 231) . and J^g is a Chern-Simons current: 

^cs = -^e^°AlF° a . (11.34) 



To derive the Chern-Simons part of flll.33j) . we used the identity d^F^ + f a bcA b ^F^ p ^ = 



and eq. ( IH.28j) . At the leading viscous order, J^g will play the role of the vorticity term in 
the hydrodynamic charge current. 

III. THE ADS BLACK BRANE 
A. The homogeneous solution 

For Anti-de-Sitter asymptotics, we choose a negative cosmological constant A = —4 in 
fill. 251) . The field equations fill. 25I) - (III. 271) then admit a family of exact solutions, describing 
a charged black brane in equilibrium. The solutions are parameterized by a set of constants: 
a 4d symmetric form h^ v with Lorentzian signature, a velocity a mass parameter m and 
a charge parameter q a in the adjoint representation of Q. In (radially shifted) Eddington- 
Finkelstein coordinates, the solution reads: 

g AB dx A dx B = -MpdxPdr + (R + rf (P^ - V(R + r, m, q a )£„Q dx^dx" , (III.l) 

A\dx A = ~^f^M^ ■ ( IIL2 ) 

The 4d indices (/i, u, . . . ) are implicitly lowered and raised with the boundary metric and 
its inverse h^ v . We may consider h^ v as a flat non-degenerate metric over the x^ 1 . The velocity 
£ M is normalized with respect to this metric, so that h^i^f = = — 1. P M „ = h^ u + £^ v 
is the projector orthogonal to with respect to h^ u . The function V(f,m,q a ) is given by: 

V(f,m,g) = l-^ + ^. (III.3) 
f = R(m, q a ) denotes the largest root of V(f, m, q), i.e.: 

m = R 4 (l + ■ (IH.4) 

This relation makes (R, q a ) more convenient than (m, q a ) as basic variables. We will find 
that these variables are related to the fluid's entropy and charge densities. 

The solution f lHLlj) - (IHL2j) is homogeneous in all but the radial direction. It does not 



depend on the gauge structure constants and couplings, or the Chern-Simons coefficients. 
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The hypersurface % defined by r = is the brane's outer event horizon. In accord with 
the notation of section III the vector £ A = (0, l^") is normal to the horizon, and we have 
9ab£ B = d A r on H. 

Substitution and direct calculation give us the horizon quantities defined in sections III Al 
EH 





— R 2 P ■ 








(III.5) 


A* 




«w = 






(III.6) 


A (0)a 


2R2 *>» 


7 -(o)m _ 

J a 






(III.7) 



Here and below, all quantities are evaluated on the horizon unless otherwise specified. The 
'(0)' superscripts anticipate the inhomogeneous corrections. To calculate 0/? and , we 
used the horizon values of the derivatives d r g^ u and 9 r A°, respectively. These are the only 
r-derivatives we will use throughout the paper; the calculation of inhomogeneous corrections 
will focus on % exclusively. In the homogeneous solution (1III.1I) - (1III.2|) . all partial derivatives 



with respect to x M vanish; also, we have f a b c A^ b A^ c = 0. Therefore: 

g(0) = Q . a W(0) =0; F^ = 0; jg^ = 0. (III.8) 
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We'd also like to have a convenient matrix of the form ( III. 61) . Let us choose A = 1/R 
and bp — so that: 

Gy,v = -^7^ - V* ■ (IH-9) 

For the solution flHI.ip - flIH.2p . this simply equals h^ u : 

G<g = V • (ni-io) 

In the decomposition (III. 1 2 j) of Q^ u , we have only the c^-term, with the value: 

cf = -kPX . (III.ll) 



B. Thermodynamic quantities 

Our black brane in equilibrium can be viewed as a thermodynamic system living in the 
4d metric h^. The velocity vector defines the system's rest frame. Using the Bekenstein- 
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Hawking correspondence between area and entropy, we define an entropy current as: 

1 



where the entropy density s is given by: 



,(0) 



-hsF 



(111.12) 



AV^h 4 

The temperature is given by the surface gravity: 

T = ^ = R_( 2 _<ht 

The energy density can now be derived as: 



(111.13) 



(111.14) 



167T 



i? 6 



3m 
16n 



(111.15) 



In equilibrium, we identify with the system's charge current J£. This will be justified 
by a conservation law in section HVBl The charge density is then given by: 



n 



An 



(111.16) 



so that: 



7(0)/* _ ,-(0)a 



(111.17) 



The chemical potential and pressure can then be derived as: 

de \ \/3<7 a 



dn a 



2R 2 



ft* A (0) I 

1 ^W\H 



p = Ts + /i a n a - e 



R 4 



q a q 



16tt V" # 6 

Finally, we define the system's stress-energy density as: 



m 
16tt 



(111.18) 
(111.19) 



-^K + (e+p)vl 



(111.20) 



Once the parameters (R, q a ) are eliminated, our knowledge about the brane's thermody- 
namics is fully described by the equation of state: 

, „s 3 /S\ 4 /3 



els, n 



7T \2 



ir 2 n n n a 



3s 2 



(111.21) 



15 



which can be derived from ( 1III.13I) . (IIII.15j) and (1III.16j) . This equation of state is conformally 
covariant, which is also evident from the tracelessness relation e = 3p. 

Using the thermodynamic variables, the horizon quantities from the previous subsection 
can be rewritten as: 



1, 



(o) 

fll' 



(4s) 2 / 3 P / 



A (0)a = _ fl a ( 



(111.22) 
(111.23) 



where c„ is again defined using the auxiliary matrix: 



G 



[IV 



(4s) 2 / 3 



1 fJLV ^fl^U ("^ 



(4 s )2/3 



with GjS = hfj, v . The other quantities can be derived directly from ( 1III.22I) 



(111.24) 



eW" = -27rT£ M r ; 



,(o) 



-hs\ 



K 



(111.25) 
(111.26) 



In the following, we will use the ansatz (1III.22|) - (1III.23|) directly, with no reference to the 
specific equation of state (IIII.21I) . or to conformal invariance. We only demand T and \x a to 
be the derivatives of some energy function e(s,n a ) with respect to s and n a , respectively. 



C. Inhomogeneous solution: outline of the calculation 

Consider again the configuration 0HI.ip - (lIH.2p . but with (h^,^, R,q a ) slowly varying 
functions of rather than constants. We will interpret this ^-dependence by treating 
(hf iV ,£ fi ,R,q a ') (or (hfj, v , i^, s, n a )) as fields on the brane's horizon H. Various quantities 
and equations can be expanded order by order in powers of the small <9 M derivatives (note, 
however, that d r derivatives are not small). We use the symbolic small parameter e to count 
these powers. We will refer to the power of e involved as the "order" of a quantity or an 
equation. 

In general, (IIII.ll) - flIII.2p with x^- dependent parameters will not be a solution of the field 
equations. However, for certain 4d configurations of (h^,^, s,n a ), it will approximate a 
solution. This has been shown for trivial, U(l) and SU(2) gauge groups, and we assume 
this to be the case for general gauge groups as well. The horizon fields of the exact solution 

1(3 



are given by functionals of (h^, s, n a ), for which eqs. (lIII.22j) - (IIII.23j) provide the zeroth- 
order terms. 

The constraints on the 4d fields (h^, i^, s, n a ) arise from the horizon projections (1II.31j) . 
( 1II.33P of the Einstein and Yang-Mills equations. We will begin by writing the first-order 
constraints for the ansatz (1III.22|) - (1III.23|) (ideal equations), followed by the second-order 
constraints using the most general corrections to f lIII.22p - flIII.23p (viscous equations). We 
require the constraint equations to take the form of conservation laws with respect to the 
metric h^ u . For the ideal equations, this requirement is automatically satisfied; however, at 
the viscous order, it restricts the possible corrections to the ansatz flIH.22p -f lHI.23p . The 
equations then describe the dynamics of a 4d fluid living in the metric h^, with configuration 
(i^, s, n a ). The null normal corresponds to the fluid's entropy velocity (since ~ 
is the area current). Most of the remaining freedom in the corrections to ( IIII.22p -( IIII.23p 
has no effect on the hydrodynamic transport coefficients. The corrections to are the 
exception, since they affect the non-abelian parts of the conductivity. Thus, we obtain a 
substantial part of the transport coefficients in a unique closed form. 

Our analysis will take place on the event horizon % of the exact inhomogeneous solution. 
In the configuration (1III.1|) - (1III.2|) with inhomogeneous (h^, i^, s, n a ), the hypersurface r = 
is null, and doesn't intersect the singularity. It is therefore an event horizon for the zeroth- 
order solution, and a zeroth-order approximation for the horizon of the corresponding exact 
solution. For the corrected solution, we will still use coordinates so that r = and g rA = l A 
on H. 

Alternatively, our calculation can be viewed from a more general perspective, as an 
analysis of the evolution equations of a null hypersurface in a hydrodynamic ansatz. The 
equations for the AdS black brane are then obtained as a special case. 



IV. IDEAL HYDRODYNAMICS 

A. Charge tensors constructed from n a 

The hydrodynamic charge density n a defines a distinguished generator in the gauge al- 
gebra. Before we proceed with the ideal hydrodynamic equations, we pause to analyze the 
simplest charge tensors that can be built from n a in a ^-covariant manner. These are given 



17 



by various products of n a with the gauge group's constant invariant tensors, such as f a bc 
and 8 a b- A reader interested only in the abelian case may skip to section IIV Bl 

1 . Scalars 

There are rank Q functionally independent charge scalars which can be constructed from 
n a . They are obtained by saturating the indices of totally-symmetric constant invariant 
^-tensors with factors of n a . We denote a full set of such independent scalars by iVj, where 
the index i runs from 1 to rankC?. The values of JVj determine n a up to gauge rotations. In 
the abelian case, the components of n a are themselves a full set of invariant scalars. 

2. Vectors 

There are rank£ linearly independent charge vectors which can be constructed from n a . 
As a basis for these vectors, we may use (iVj) a = dNi/dn a . Because they are constructed 
from n a and constant invariant tensors, the (Nj) a all commute with n a and with each other 
under the gauge algebra. For non-singular values of n a , the (Ni) a form a basis for the 
centralizer Z(n) of n a , i.e. the subspace of the gauge algebra which commutes with n a . 

3. Matrices 

One type of charge matrix which can be built out of n a is the external product (Ni) a (Nj) b 
of two vectors. For generic values of n a , such matrices span the entire matrix space over 
Z{p). In particular, their linear combination gives the projector Z a b onto Z{n) with respect 
to the Killing metric (for an abelian gauge group we define Z a & = <J a &, and similarly for 
abelian pieces of a semi-simple group). 

There are also charge matrices which do not lie in Zip), and are not reducible into 
products of covariantly constructed vectors. To clarify the structure of such matrices M a b, 
note that for any vector (Ni) b , the expression M a b(Ni) b is a charge vector constructed from 
n a . Therefore, M a b(Ni) b lies in Zip). It follows that for generic values of n a , the matrix M a b 
can be expressed in block- diagonal form: 

M ab = Mi + Mi, (IV. 1) 
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where both indices of M ab and M^ b lie respectively in Z(n) and Z ± (n), the orthogonal 
complement of Z(n). The projector Z^ b = 5 ab — Z ab onto Z^ip) is an example of a M^-type 
matrix. Another trivial example is / a b c n c . 

The matrix M ab is constructed only from n a and invariant constant tensors; therefore, it 
must be invariant under rotations generated by n a . The corresponding infinitesimal state- 
ment is: 

fa dc n c M db + f bdc n c M ad = f adc n c M db - M ad f dbc n c = . (IV. 2) 

In other words, M ab commutes with f abc n c in the matrix sense. For the M^-piece, this is 
trivial. For the M^-piece, the commutation implies that it must be a matrix polynomial in 
fabcn c - This is because for generic n a , the matrix f a bcn c has distinct eigenvalues over Z ± (n), 
since are no multiplicities in the nontrivial roots of compact simple Lie groups. But if any 
M^-type matrix is a polynomial in f abc n c , then these matrices all commute with each other. 
Since they also commute with all M^-type matrices, we have: 

Proposition IV. 1. All matrices constructed from n a with indices in 2'" L (n) commute with 
all other matrices constructed from n a 

B. Current conservation 

At first order, the restricted Yang-Mills equation ( III. 331) can be written as: 



+ f abc {Af b j^ - j (0) ^4 1)c ) = 0(e 2 ), (IV.3) 

where "(1)" superscripts indicate first-order corrections to the zeroth-order expressions 
(IITL231 . 



1. Abelian case 

For abelian Q, the / a f, c -terms in (1IV.3I) vanish, and we immediately get the ideal current 
conservation equation: 

= 0(e 2 ). (IV.4) 
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We can now substitute the partial derivative with D^, the covariant derivative with respect 
to h^y. Plugging in the zeroth-order value from section flHIBI) and dividing by y/—h, the 
equation becomes: 

Dn a + n a D^ = 0{e 2 ), (IV.5) 

where D = l^D^ is the directional derivative along i^, and = P^ u D^i u is the fluid's 

expansion rate. Eq. (1IV.4|) justifies our interpretation (1III.17I) of jl ^ as the zeroth-order 
conserved charge current Ja M . 

2. Non-abelian case 

For a non-abelian charge group, the same ideal conservation law is obtained. First, let 
us expand the / a 6 C -terms in (1IV.3I) as: 



/a 6c (4 0)6 i (1)cM - jW»Ajp*) = -f abc (»%j (1)cti + n b FA^ c ). (IV.6) 

Let us contract this with a charge vector (Ni) a which commutes with n a . The second term 
in (IIV.6j) vanishes immediately. As for the first term, note that \x a is a (/-covariant function 



of n a and the scalar s. It is therefore in the centralizer Z(n) of n a , and commutes with 
(Ni) a . We conclude that contraction with (iVj) a annihilates both terms in f ]IV.6j) . Since the 



vectors (Ni) a span Z(n), we may now write the Z(n)-projection of (1IV.3I) as: 

Z a %j^ = 0{e 2 ) (IV. 7) 

Z ah Dn b + n a D^ = 0{e 2 ). (IV.8) 

Now consider the Z J -(n)-projection of ( 1IV.3I) . It describes the vanishing of a first-order 
geometric scalar with a (dimQ — rank component charge index in Z J -(n). Now, the only 
independent quantity of this type which can be constructed from (h^,^, s,n a ) is Z^ b Dn b . 
Thus, unless the corrections A^ a and jl 1 ^ are such that eq. (IIV.3j) becomes degenerate, it 
necessarily implies: 

Z^ b Dn b = 0(e 2 ). (IV.9) 

Conversely, suppose that non-degenerate first-order hydrodynamic equations exist (whether 
or not they are fully captured by the horizon equations flll.31 j) . flll.33j) ). The preceding 
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discussion shows that these equations must contain eq. (1IV.9I) . But then, that relation itself 
implies that the / a b c -terms in flIV.31) must vanish. Indeed, f a bcH b £ f ij^ CfJ ' and f a bcn b ^A^ c are 
first-order scalars in Z J -(n); if eq. ( IIV.9P holds, then these quantities vanish at first order. 

We conclude that the projected Yang-Mills equation (IIV.3P takes the form of the conser- 
vation law (jrV.4j) - flIV.5j) in the non-abelian well. 



C. Gauss-Codazzi equations and energy- momentum conservation 

The Gauss-Codazzi equation f ]II.31j) reads at first order: 



c (o) duS (o)u + 2S^>d [v c^ = l^j^d^f + 0(e 2 ). (IV.10) 

Plugging in the zeroth-order values from section IIII Bl and substituting <9 M with D^, the 
equation becomes: 

-^^h(D u {Tst^t) + sD^T) 

(IV.ll) 

= (D^ny^) + n a D^ a - iiJ^D v (n a t)) + 0(e 2 ). 

In the last term on the RHS, we recognize the current divergence d u j^ u . Using the ther- 
modynamic identities e + p = Ts + fi a n a and dp = sdT + n a dfi a , we rearrange the equation 
to get: 



-8kD v T<?> = -Sir^iM^ + 0(e 2 ), (IV.12) 



where Tji ^ is the stress-energy density defined in (1III.20|) . Using eq. (1IV.4|) and dividing 



by — 87r, we get an equation of hydrodynamic stress-energy conservation: 

D u T^ )u = 0{e 2 ). (IV. 13) 

Together with the current conservation law (1IV.4I) . we now have a complete set of equations 
for ideal charged hydrodynamics. Note that these equations are independent of the structure 
constants f a f, c and the Chern-Simons coefficients f3 a b c ; these represent interactions between 
adjacent regions of the brane/fluid, which will come into play only at the next order. 
In components, eq. (1IV.13I) reads: 

De+(e + p)D IM e = 0(e 2 ), (IV.14) 
(e + p)Dt + P^D v p = 0{e 2 ). (IV.15) 
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A convenient rewriting of eqs. (1IV.5I) and (1IV".14j) is: 



Ds = -sD^ + 0(e 2 ), (IV.16) 

71 

D— = 0(e 2 ). (IV.17) 



s 

Eq. flIV.16j) expresses the conservation of entropy in the evolution of an ideal fluid: 



= 0{e 2 ). (IV.18) 
It can be derived more directly from the focusing equation (III. 321) . which reads at first order: 

k6 = 0{e 2 ) -d^ = 9 = 0{e 2 ). (IV.19) 

This tells us that in the ideal approximation, the brane evolves in an area-preserving manner. 



V. FIRST-ORDER CORRECTIONS AND VISCOUS HYDRODYNAMICS 
A. Classification of first-order quantities 

We pause to consider the possible first-order quantities which can be constructed out 
of (h^ u , s, n a ). Some quantities are first-order for generic configurations, but become 
restricted to 0(e 2 ) by the ideal constraint equations; here we discuss quantities which remain 
first-order after the ideal equations are imposed. 

1 . Scalars 

Eqs. (liV.16j) - flIV.17j) imply that is the only independent first-order geometric scalar. 
Other commonly encountered geometric scalars have the form Df, where / is some function 
of state (which may carry unspecified charge indices). Using flIV.16j) -f )IV.17j) . such quantities 
can be related to as: 

D f = -*(%) j D,t + 0{e 2 ). (V.l) 

In particular, all first-order scalars with a single charge index have the form f a D [li l li , where 
f a is a charge vector constructed from (s, n a ). As a result, the charge index of such quantities 
necessarily lies in Z(n). 
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2. Vectors 



Let us list the possible first-order vectors transverse to i^. First, we have the acceleration 
Di^ and the transverse gradients P^ u d u f of thermodynamic functions / (with possible charge 
indices). This list is redundant, since eq. (1IV.15P relates D^ 1 to P^ v d v p. Also, eq. (1IV.17P 



implies that d fl (n a /s) is already transverse, with no need for a projector. Thus, a general 
transverse gradient can be written as: 

Once axial quantities are allowed, we also have the vorticity vector density: 

^ _ ^u P „ LdpL = l e ^ iuUJpa 9 (V.3) 



where is the vorticity tensor defined in (IV".6[) . Finally, we can get a first-order vector 



along by multiplying it with a first-order scalar. 

We conclude that the most general first-order vector has the form: 

y{i)n = vfPDJT + VD^s + V a D»— + ^=u» , (V.4) 

s \/—h 

where the coefficients V, V, V a and V are (/-covariant functions of (s,n a ). To get a charged 
first-order vector, one should add charge indices to these coefficients in the obvious manner. 

3. Rank-2 tensors 

The only first-order traceless rank-2 tensors transverse to are the shear tensor ir^ and 
the vorticity tensor u^ iu : 

1 

3 



V, = P^D {p e a) - -P^D/o , (V.5) 



V — - 1 fx 



P^D [p £ a] = -Un [p d/ U] . (V.6) 



Other first-order rank-2 tensors can be constructed from the tensor product of with 



first-order vectors, or of P pv with first-order scalars. 



B. Corrections to the horizon ansatz 



Before we write down the second-order constraint equations, we must consider the pos- 
sible first-order corrections to the horizon ansatz (lIII.22l) - flIII.23j) . These corrections can be 
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restricted by a first-order fixing of variables. We consider s,n a ) as the basic variables 
in (IIII.22I) - (IIII.23I) . with T and \i a derived from (s,n a ) through the equation of state. We 
fix the direction of the entropy velocity £ M and the magnitude of the entropy density s by 
requiring that eq. (IIII.25P remains without corrections, i.e.: 

= o ; „(!) = o . (V.7) 

Due to this condition, the correction to 7 M „ must be transverse to I 11 and traceless with 
respect to h^ y . We will find that the precise form of 7^ is otherwise irrelevant to the 
constraint equations. For completeness, we note that this correction must take the form: 

7$ = , (V.8) 

with a some function of (s,n a ). 

Let us now turn to the 7 MV -independent components c M of the extrinsic curvature. For 
the definition of c M , we still use the auxiliary matrix (jTTL24j). The most general correction 
to c M reads: 

C W = cl^DX + cD^s + CaD^— + -S=c^ , (V.9) 

s v — n 

where c, c, c a and d are some functions of (s,n a ). The requirement for the constraint 
equations to take the form of conservation laws will place restrictions on these functions. 
Otherwise, they will have no impact on the hydrodynamic transport coefficients. 

We fix the meaning of the charge density n a by requiring the relation i^j^ = —\J—hn a 
to hold without corrections: 

«^ = 0. (V.10) 

As we will see, this condition fixes n a as the longitudinal component —l^J^jyJ—h of the 
conserved hydrodynamic charge current. As implied by ( III. 2411 . all the remaining first-order 
corrections to j% are determined by the horizon gauge field A a . The general correction to 
Af, takes the form: 

A# = AJ^D u t + A a D,s + A ab D^- + , (V.ll) 

s \J—h 

where A a , A a , A a b and A' a are some functions of (s,n a ). Again, the requirement for a 
conservation-law form of the equations places restrictions on the form of these functions. In 
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the abelian case, they otherwise have no impact on the resulting transport coefficients. In the 
non-abelian case, however, the Z J -(n)-components of A a b play a crucial role in determining 
the Z J -(n)-piece of the conductivity matrix. 

C. Current conservation 

In this subsection, we analyze the second-order terms in the projected Yang-Mills equation 
(III. 331) . The abelian case is simpler, and we will consider it first. 

1. Abelian case 

When f a b c = 0, the horizon's intrinsic field strength F^ u can be written as: 

F; v = 2d { ,A^ a + 0(e 2 ) = 2D lv (ji't A ) + 0{e 2 ). (V.12) 
To find ja , we write down the first-order terms of eq. (III. 24ft : 



l^j {0)au + lf}j {l)av = . (V.13) 

The first term vanishes due to 7^L* = 0. The other terms evaluate as: 

(4s) 2/3 P^ji 1)U = {P;D^ a + HaDiJ. (V.14) 

7T ^ 



Together with eq. (IV.lOj) . this gives: 



iP" = - ^t^ KD^a + VaD£»). (V.15) 

From equation (1IV.15I) and the identity dp = sdT + n a dfi a , we derive an identity of ideal 
hydrodynamics: 

P;D^ a + ^Dl, — T (s h a — t^L\ P;D V ^ + 0(e\ (V.16) 
which enables us to rewrite ( IV.15P as: 



,6 



= -^hTa ab P;D u ^- , (V.17) 



T 

'>b_ 

2 4 /% V ao e + P 



& ab = T^FTZ Sab - ) • (V.ll 
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For the Chern-Simons current (111.340 . we have: 

jg)a„ = _^abc e ^ P a A (0)b d A (0)c = _ £ ^ 1 V ,y lg) 

71 7T 

To second order accuracy, the horizon Yang-Mills equation ( 1II.33j) now reads: 

<V£ = 0(e 3 ), (V.20) 
where the conserved current is given by: 

= V^h (n a F - Ta ab P» v D u ^) - - f3 abc + 0{e 2 ). (V.21) 
2. Non-abelian case 

The non-abelian case requires a more elaborate argument. The first-order field strength 
F$ a now receives an additional contribution: 

2f\ c A^ b A^ c = -2f a bcf i%A^ c . (V.22) 

The factor of f a bc^ b singles out the Z J -(n)-components of Al 1 . Of all the terms in ( IV. 11 II . 
this leaves only the v4 a b-term; the other coefficients are charge vectors built out of (s,n a ), 
and therefore lie in Z(n). Furthermore, we may rewrite the A a t,-term as: 

n b 

A ab D^- = A ab D^ b + ... (V.23) 

where the ellipses indicate terms consisting of (zeroth-order) charge vectors in Z(n) mul- 
tiplied by (first-order) charge scalars. These terms will again be annihilated by f a b c ^ b - In 
summary, we may write: 

iV = 2 (ZatD^e^) + j^D v]f i b ) + 0(5 2 ), (V.25) 

where is the following matrix, whose indices lie in Z J -(n): 

ji = Zi - f acd ^A\ . (V.26) 
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Analogously to (IV. 17ft in the abelian case, we can now find the first-order correction to j£ 



as: 

,1/3 



A 1 * = ~^hT (ai + ^-j^j P<» Dv (L , (V.27) 

where we used j a L b /i b = to manipulate the j^-term, and the matrix a a \ with indices in Z(n) 
is defined similarly to ( 1V.18j) : 

^ = 7??) ■ < v ' 28 > 

Turning to the Chern-Simons current (III34p . it's easy to see that the first-order expression 
dVTT9|) is valid in the non-abelian well: there is no contribution from / a (, c -terms, 



because A^ )a ~ ip. Examining eq. rtVTlQl) . we see that (3 a bcf^ b ^ c is a gauge vector covariantly 
constructed out of (n a ,s). It follows that (3 a b c fi b fi c , and therefore J^g itself, lie in Z(n). 

We are now ready to write down the second-order terms in the LHS of the restricted 
Yang- Mills equation fill. 33 j) : 



7T MP " 



7T 



n 



(V.29) 



Let us contract (1V.29|) with a charge vector (iVj) a covariantly constructed from n a (recall 



that such vectors span Z(n)). We find that all the terms in the last two lines vanish: 

• In the j£g CAt -term, we have f a bc(Ni) a Jcs Cfl = 0, because (Ni) a , J^s^ e Z(n). 

• In the (j( 2 ^ + 4 2 s )cM )-term, we have f a bc ( N i)aH b = 0, because {Ni) a ,fi b e Z(n). 

• In the -term, we have f a bc(Ni) a n c = 0, because (iVj) a G Z(n). 

• In the /3 a f, c -term, we have f bde fid P a bc(Ni) a iJL c = 0, because (3 a b c (Ni) a n c E Z(n), as a 
charge vector covariantly constructed from (n a , s). 
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In the contraction of the f a bc A,j^ CIJ, -teTm with (Ni) a , only the Z J -(n)-components of A^ b 
and survive. We get: 

,/ — h « 1 / 3 T n e 

(N t ) a r bc A^ b j^ = - 2 , /H r b \N t ) a A bdJ iP^D^ d D^ . (V.30) 
Generically, f a bc(Ni) c is an invertible matrix over Z (n). We may therefore denote: 

Jab = K'fcbdW = fac d (N t ) d k C b , (V.31) 

where we used the matrix commutation of k ab and f a bc{Ni) c (Proposition IIV. 1]) . Expression 
f ]V.30j) now becomes: 

J —h s 1 / 3 T - n e 

-f ab %NA a A hd kJf feg {NA 9 P^D v ^ d D,^- 



(V.32) 



24/3^ J \ ' / / cc J oct c jjeg\+<i) * 

= f V 24/% r bc {N % ) a A u kJ j fe9l fP^D vl i A 

/ A /Z/7 s i/3 _ \ 

= - (iV;)^ (^ V 24/37r Aujifffeg^P^D^j 

= WD, { ff^ T {jt - JcJ^D^ • 
Thus, the full contraction of (1V.29|) with (iVj) a reads: 

(-v^T (4, + ^(2^ + jiji)) P""D W £ + J«^) • (V.33) 

The j^-term is actually redundant. Indeed, let X ab = X[ ab ] be an antisymmetric matrix with 
indices in Z J_ (n). As our spanning set (Ni) a of Z(n), we may choose the derivatives dNi/dfi a 
of a set of scalar functions of /j a . Then (Ni) ab = d(Ni) a /djj, b is symmetric^. Furthermore, 
according to Proposition lIV.lt X a b and (Ni) ab commute as matrices, so their matrix product 
is antisymmetric. We then find: 

WD, (x ab P^D u ^) = WD, (^X ah P^D v A = -^X ab P^ u D,(Ni) a D u /j, b 



(V.34) 



rj-\ j \ <> / \ J 

~X ab (N i ) a c P^D, f i c D„fi b = 



2 Since the (Ni) a span Z(n), this symmetry property follows for every charge vector constructed from /x a 
(or n a ). 
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where we used the relations X a bfi b = X a b(Ni) b = and the antisymmetry of X a b(Ni)^. We 
conclude that the j^-term in flV.33j) may be replaced with any other antisymmetric matrix 



over Z J -(n). We can then write (1V.33P as: 

{Nd a D, (-V=hT(*i + atiP^D^ + jg}<") , (V.35) 
where a^ b is a function of (s,n a ) with indices in Z ± (q) whose symmetric piece is given by 



s 



1/3 



a (ab) — l^/5Z3ca3d> ) (V.36) 

and the antisymmetric piece is otherwise arbitrary. Eq. flV.36j) implies that ah^ is 
positive semi-definite. As we will see, this property is related to the second law of thermo- 
dynamics. 

Since the (iVj) a span Z(n), we conclude from ( ]V.35j) that the Z(n)-projection of the 
second-order terms (lV.29j) is the Z(n)-projection of a divergence. Combining this with the 
results of the previous order, we get: 

Z ab d^ = Z ab d^(4 0)fl + J^) = 0(e 3 ), (V.37) 

where Ja ^ is given by 

= _ T v ab p^ Dv t- _ _/3 a6c/i y^ + 0(e 2 ) (V.38) 

1 IX 

with a ab = al b + a^ b . 

Let us now return to the full second-order expression (1V.29|) . We see that it depends on 
second-order corrections to the horizon ansatz. Specifically, this dependence is contained 
in the terms f a bc(A^ b j^ Cfl + A^ b j^ c>1 ). At first sight, we may be concerned that this 
quantity is sensitive to the second-order fixing of variables. Reassuringly, this isn't so. 
Indeed, consider a second-order redefinition of the hydrodynamic variables (h^,^, s,n a ). 
A redefinition of variables means that the value of geometric and gauge quantities such 
as f a bc A b ^ must remain unchanged, while their functional dependence on (h^, i^, s, n a ) 
may be altered. In particular, the contribution to f a bc(A^ b j^ Cfl + A^j^ 01 *) due to the 
redefinition must equal minus the change in the value of f a b c A^ b j^ c>l = \J —hf a b c /j, b n c . 
But f a bc^ b n c always vanishes regardless of the variable-fixing, since /i a is always a covariant 
function of (s,n a ). We conclude that expression (IV. 291) is insensitive to the second-order 
fixing of variables. 
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Now, since the matrices f a b c n c and f a bc^ b are invertible over Z (n), we see that appro- 
priate choices of A^ a and jjj 2 can make the Z J -(n)-projection of (IV. 291) equal anything. 
In particular, we can make it equal Z^d^J^ 1 ^ , with any choice of <r^ in (IV.38I) . Then eq. 
(IV.37P is upgraded into a full conservation law: 

<V£ = 0(e 3 ). (V.39) 

Thus, for an arbitrary first-order correction A^ a to the gauge potential, there exist ap- 
propriate second-order corrections (Af, ; , ja ) for which the horizon Yang-Mills equation 
takes the form of a conservation law. The first-order correction determines the symmetric 
part of the Z J -(n)-piece of the conductivity matrix, while the second-order corrections 
determine its antisymmetric part cr^u- 

One of the allowed sets of first-order and second-order corrections corresponds to the non- 
abelian AdS black brane with vanishing gauge fields at r — > oo. Unfortunately, we cannot 
find the relevant corrections on the horizon without solving the radial equations. Thus, 
we do not obtain the specific form of a^ b . However, our analysis provides two non-trivial 
relations between the gauge fields on the horizon and on the AdS boundary: 



a (ab) a PP ears m the conserved current (TV38j) . which should correspond to the boundary 



quantity ( 11.61) . On the other hand, ah^ is defined in ( 1V.36I) in terms of j^ b , which is 
derived from the gauge potential A^ on the horizon. This provides a first-order 
relation between horizon and boundary fields. In the simplest non-abelian case Q 



SU(2), we can check this relation against the results of the bulk calculation in 26]. 



In the notation of 



the first-order relation we find reads: 



lim r2 9 W = ll (-2 + 2 / (1) + ? 2 ((/ (1) ) 2 + ^ (1) ) 2 )) 



(V.40) 



H 



where f^ 1 ' and g^ are the coefficients of the two non-abelian terms in A", and the RHS 
is evaluated at the horizon. This is a scalar equation, because for an SU(2) group, 07^ 
has a single component ~ (q 2 5 a t, — q a Qb)- Following the numerical prescription of [26 1 
for solving the radial differential equations, we have calculated f^(r) and g^'(r) for 
several values of q 2 /R 6 in the range 0.01 — 1.7, and found that (1V.40I) holds within the 
numerical accuracy of ~ 1%. We view this as evidence for the correct identification of 
the conserved current (IV.38j) with the AdS/CFT boundary current (II. 6j) . 
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• The condition that the Z ± (q)-projection of (lV.29j) should equal Z^d^J^ 1 ^ provides 
another relation, at second order, between horizon fields (A^ a , A^ a ) and boundary 
fields (J^ 1 ^, with its and cr^ terms). 

D. Gauss-Codazzi equations and energy-momentum conservation 

In this subsection, we analyze the second-order Gauss-Codazzi equation. The non-abelian 
case doesn't lead to significant complications, so we will handle it straight away. 

The horizon's first-order shear/expansion tensor can be derived directly from the zeroth- 
order metric via ( 111.21) : 

V = Unj?J + 0(e 2 ) = \£»D pl $ + 7 $Ao* p + 0(e 2 ) 

= {As) 2 " U, u + Ys [Ds + sD ^ P A + °^ = ( 4s ) 2/3 ^ + 0(e 2 ). 

In the second equality, we wrote the Lie derivative in terms of the connection D^. In the 
last, we used the ideal equation ( 1IV.16I) . Decomposing the result with ( 1II.4p -( IIL5l) . we get: 



0(1) = 0; = (4 S ) 2 / V ■ (V.42) 

The result for was already known from the ideal equations. 

Let us now evaluate the second-order terms in the Gauss-Codazzi equation fill. 3 1 j) . For 
the LHS of flll.31j) we get, using our choice A = l/(4s) 2//3 : 



((^e^r) + 2^cg ) , (V.43) 

where we can use the zeroth-order expression (IIILICJ for G^ u , substituting (G~ 1 ) pv with h pu 
and with D u . This gives: 

Dv ( J^ e p V ) + = D » ( 4 V^K) + 8v^r^c« . (V.44) 
We now turn to the RHS of (111.3ip . At second order, it reads: 

to{F^jf> v + Fj&£>). (V.45) 
The first term can be expanded as: 



167iV^hn a rd [fl A^ a + 0(e 3 ), 



(V.46) 
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where we used the fact that both A ( p ' a and the first-order geometric scalar £ v Ah) lie in Z(n). 
The second term in (1V.45I) reads: 

toF$*fi> = -Svrv^ (t^P» p D p ^ (!>„(/,%) - fD^Q 

(V.47) 



where we used eqs. (fV25|) . (TV27j) . ( TV6]) and (TVT36]) . and the fact that af ab) ^ b = 0. 

There is another second-order correction that must be taken into account. At first order, 
we brought the Einstein-Gauss-Codazzi equation to the form (HV.12j) . The RHS of that 
equation vanished due to the first-order current conservation. At second order, this is no 
longer true. Instead, we have: 



-8^ a £,D v (r a -J^) + 0(e s ) 

Sir^DvJ" - 8n£ p (n a D v ( ' V^hTa l } b P up D p ^ + ^ ahc ^ b ^ c u v 



T 7T 



-haiP up D vt i a D p fi b )+0{e 



(V.48) 



where we used eq. ( 1V.38I) and applied the Leibnitz rule on the cx^-term. We see that the 
cr^-terms in (IV.47|) and (jV.48p cancel each other. The cxjj b -terms combine to give: 



8nD„ V-hT f , a a^P vp D i 



' T 



8nV-hT^ b D^D,£ u . 



V ^D,£ V . (V.49) 



We must now use the Leibnitz rule repeatedly to manipulate the second term of (IV. 49j) . 
keeping the products Tfi a a^ b and /i fe /T intact. The expression becomes: 



^hT^J [p P p v) D p ^ 



+ Svrv^ (d v (p;T^D^ - D^D p {T^ ab ) + D u {T^l b t)D 



V rp 



(V.50) 



The first line is the divergence of a symmetric tensor; it can be used as part of a stress-energy 
conservation law. The terms in the second line we intend to cancel with other contributions. 
For that purpose, we decompose the gradients in the last two terms into combinations of 
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d^s and d IJi (ji a /s), and use the ideal equation (IV. ip . The second line then reads: 



V \ , m . a-J\ n „u fd(V b /T) 



(V.51) 



dn c 7 C T \ dn 




We now turn to the Chern-Simons term in (IV.48I) . To approach it, we will need two 
identities concerning the vorticity: 

= 2u»Dl p , (V.52) 

u"u) uil = . (V.53) 

The first of these is easy to derive from the definition (IV.3P and the four- dimensionality of 
the horizon. To derive the second, note that u pu is effectively three-dimensional, as it is 
transverse to This implies that w^ui p(T \ = identically. Eq. ( TV53|) then follows from 
the definition (IV. 31) . Written out more explicitly, (IV. 53|) reads: 



u v (DJ„ - D p l v - £ p DQ = . (V.54) 

Using these identities, the index symmetry of (5( a bc) and the ideal equation (IIV. 16[) . we can 
expand the Chern-Simons term in (IV.48I) as: 

64 n „,o a b c„ ( /WV/^ (V ' 55) 

The first term will become part of the stress-energy conservation, while the second will be 
canceled by other contributions. 

It remains to analyze the Z u d [u c { ^ term in ( rV44|) and the td [u A { ^ a term in dV46|) . They 
will be used to cancel the non- divergence pieces we picked up in ( IV.5ip and (lV.55p . First, we 
gather both terms on the same side of the Gauss-Codazzi equation by moving the td [u A^ a 
term to the LHS. We then plug in the expressions (IV.9P and (TV lip for qp and Aj^ , and 



33 



use the ideal equation (IV. lj) to get: 



V^h{std [v c^ + 2im a td [u A ( $ a ) 

= A^h ]D V (p;(sc + 2nn a A a )D p £^ - (c + 2nn a A a /s)D v £ u D p s 

+ . h (*±|^M) a - s ( a(B - + y jW ' ) ) - ma. + 

+ 8sfL> [i , (w H (c / + 2v™X/s)) • 

(V.56) 



This expression cancels with the non- divergence terms in (IV.51j) and (IV. 55j) if the coefficients 
in (jV3|) and dVTT]) satisfy: 

sc' + 27in a A' a = ^ Qbc /i>V c (V.57) 
sc + 2nn a A a = -2nsT/i a a l } b ( ^J^ ^j (V.58) 



n/s 



(9n a 



en. ( 



(V.59) 



These constraints on (c, c, c , c', A a , A a , A a b, A' a ) clearly have non-unique solutions. Further- 
more, they place no restriction on the Z J -(n)-part of A^, which affects the symmetric 
Z ± (n)-piece of the conductivity matrix. 

Finally, we should ask whether the £ v d[ v cty and £ u dy p A^ a terms can not only cancel 
unwanted pieces as described above, but also generate additional contributions with the 
form of a symmetric tensor's divergence. If so, such contributions could be included in the 
hydrodynamic stress-energy conservation. The answer, however, turns out to be negative. 
On one hand, the £ v d^ v c^) and £ u d[ u A^ a terms are necessarily transverse to £^. On the 
other hand, consider the most general first-order symmetric tensor density: 



V^h (riTT^ + D p £ p (r 2 P^ + r 3 ff) + 2t a &D£ v) + 2t£&D u) ^ 



(V.60) 
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The divergence of such a quantity is found, using the ideal equations, to have the following 
component along l^: 

- (r 3 + r A )P^D^Dt u - Dip ((r 3 + 2r 4 )D^ + D"r 4 ) - r 3 i2<^r (V.61) 

77 77, \ 

- D^(D^ + 2r 5 W) - r^D,D v ^ J - w"(D„t 6 + 3r 6 D^), 

where i?^? is the Ricci tensor associated with the metric h^. From considering the indepen- 
dent terms in ( IV. 6 111 , we see that the whole expression can vanish only if all the coefficients 
(ti, . . . , Tq) vanish. We conclude that there is no symmetric first-order tensor whose diver- 
gence is transverse to -P. Thus, one cannot generate the divergence of such a tensor from 
t v d[ v (y~) and fd[ v A*~) a by choosing appropriate first-order corrections. 

Collecting our results, we have brought the Einstein-Gauss-Codazzi equation to the form: 

-&>*D v {Tp> + T^>) = -8n^D„r a + 0(e 3 ), (V.62) 

where the first-order stress-energy density is given by: 

T (Da- = + 2Tfi h al&P^D p ^) - l^^VV^^"' 

\Zn 1 J 671 

= -V^h ( -Lit" + ^ - &P^dM (V.63) 

Using the current conservation ()V".37j) and dividing by — 87r, eq. ( 1V.62j) becomes an energy- 
momentum conservation equation: 

D v (Tp> + T^>) = 0(e 3 ). (V.64) 

The viscous stress-energy density T^ u + T^ u can be written as: 

= (phT + (e + p)u»u u - — tt^) + 0(e 2 ), (V.65) 
V 2n / 

where we defined the energy velocity = + 0(e) as the timelike unit eigenvector of T^ v . 
Explicitly, u M is given by: 

e + p T 3n(e + p) 

sl/3T ( ..a W b n a \ nuUT ^ fj, a 4f3 abc fi a fi b fi c u ( ^ r 2 . 



2 4 / 3 7r(e + p)V e + P / T 3vr(e + p) v ; 
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At the given order, it makes no difference whether ir^ and are defined in terms of £ M or 
in terms of u M '. Using rather than as the basic velocity variable produces the equations 
in the Landau frame. We use (1V.38|) and (1V.66|) to write the conserved charge current in 
terms of 

where the conductivity matrix is: 

°ab = (Tab + ^ab , (V.68) 

(V.69) 



|| _ -|| n a ^ c _\\ s 1/3 / /i c n a \ / ' [x c n h 

a - b = a " b ~ 7Ti* ~ ¥Jh \ Zca ~7T~p)\ b ~7T~ P 



sl/3 ( rr ^( a n b ) fi c /J, c n a n b 

^ab ! r 



24/3^ e + p (e + p) 2 

The relation ( 1V.66|) between and can be rewritten using a\ h : 



s T 37r(e + p) 



Finally, from this we derive the expression for the entropy current in the Landau frame: 

,a ,,b , ,c 



(V.71) 



E. The focusing equation and dissipative entropy production 

It is a standard exercise to derive the second-order entropy production rate d^s^ from 
the viscous hydrodynamic equations. Alternatively, it is convenient to derive entropy-related 
results from the focusing equation ( III. 3211 . which deals directly with the rate of area produc- 
tion on the horizon. We will now demonstrate this approach. For the RHS of ( III. 3211 . we 
have, using (TVT47L (TvUBD . (jV\28j) , and ( IVT68|) - fl\T69|) : 

^^j: = ^^^^ + o(. 3 ) 

_ 2vrT 2 / S V3 / Ha n b) fi c ^n a n b \ ,\ f ^ b 3 

- — \¥i^\ Zab -^ + jrTW) ab ) D »T D »T + ° i£) 



ItxT 2 u a u b 

Z ' H1 _ nui/n t—ni— 

' jl rp rp 



-a ab P^D^D u ^- + 0(e 3 ). 



(V.72) 
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At second order, the LHS of (III. 321) reads: 

K (0) - — ^C^V^A" = ^ Td - + ( V - 73 ) 

where we used our choices G$ = h^ v and A = l/(4s) 2 / 3 , and substituted (IV.42|) for affi- 
The focusing equation then becomes: 



2nT6 - tt^ = a ah P^D^D u ^- + 0(e 3 ) (V.74) 



2ttT 2 u a u b 

V" — ~ a abP D^ — D v — 

=► 9 = la ab P^D^D u ^ + 2^V^" + 0(s 3 ). (V.75) 
Both pieces of <j( a b) = cr^ b + a ua>) are positive semi-definite, as can be seen from ( IV. 6911 and 



(]V.36|) . As a result, the RHS of (jV.75j) is non-negative. The rate of area production can be 



derived from ( TV.75|) as: 

d ^ = v9 = «<°>0W + 0(e 3 ) = V^h U T a ab P^D^D v ^ + ^W^) + 0(e* 

(V.76) 

which translates immediately into an entropy production rate: 

V = \d»S» = V^h ( T a ab P^D^D„^ + ^ w"") + 0(8'). (V.77) 

VI. DISCUSSION 

In this paper, we've analyzed the dynamics of the event horizon of a boosted Einstein- 
Yang-Mills black brane. The corresponding equations define the relativistic viscous hydro- 
dynamics of conformal field theories with non-abelian conserved currents. We introduced a 
non-abelian Chern-Simons term and derived from the horizon dynamics the hydrodynamic 
constitutive relations in the presence of anomalous non-abelian global symmetries. The null 
nature of the event horizon plays a crucial role in the success of the derivation: the amount of 
independent 4d projections of the bulk fields is reduced by the relations (III. 101) and (III. 241) . 
which are unique to null hypersurfaces. 

As we have seen, the calculation lends itself to a more general context: it may be applied 
to the dynamics of various null hypersurfaces with Einstein- Yang-Mills fields, in an ansatz 
where the hypersurface's evolution equations take the form of hydrodynamic conservation 
laws. Equations of state other than (1111.211) can be used. Presumably, different equations 
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of state correspond to different equilibrium bulk solutions, with eq. ( 1III.21I) arising in the 
special case of the AdS black brane. Once an equation of state is chosen, in the non-abelian 
case there is still the freedom to choose the -Z' ± (n)-piece of the conductivity a^ b . The value 
of cr^b arises from the corrections to the bulk gauge potential A a ^. The standard choice in 
AdS/CFT arises from the condition (II. 5p . Given an equilibrium bulk solution, different 
inhomogeneous corrections to A ^ lead to different non-abelian conductivities for the same 
equation of state (provided the corrections are such that the projected Yang-Mills equation 
can be cast conservation law). 

On the other hand, we see that Einstein- Yang-Mills fields on a null horizon cannot encode 
an arbitrary hydrodynamic system: once the choices discussed above are made, the other 
transport coefficients are uniquely determined. Specifically, the shear viscosity, the bulk 
viscosity and the Z(n)-piece of the conductivity matrix are fixed by the null hypersurface 
equations to specific functions of state, while for a general fluid they may be arbitrary (under 
the restriction of positive semi-definiteness). 

In addition, we obtained a unique form for the vorticity coefficient in the current (1V.38j) . 



According to the argument presented in 



22], this coefficient is not arbitrary for a general 



fluid, but uniquely determined from the anomaly coefficients of the underlying; field theory. 



Our calculation produced the exact form for the coefficient prescribed in 



22|, with anomaly 



coefficien ts g iven by C a bc = — (Z/^Pabc- Furthermore, our result directly generalizes the 



result of 



22( to the case of non-abelian charges. We intend to expand on this subject in a 



separate work. For now, we briefly note that the argument in [22] can be carried through for 
non-abelian charges by replacing all spacetime derivatives with gauge- covariant derivatives 
with respect to the external fields. 

Our non-abelian results reduce to the abelian case in two separate ways. First, we have 
the case where the entire charge group is abelian to begin with. Then all charges commute 
with each other, and Z L {ja) is the zero subspace. We therefore substitute Z a b = 5 a b and 
a ab = in our formulas, reproducing the known results for the abelian AdS black brane. The 
abelian limit should also be obtained for small charges and for weak couplings. Neglecting 
the /abc-terms, eq. (IV.26j) becomes = Z^ h , while in flV.29j) only d tl (j^ a>l + Jqq) remains. 



Plugging in eq. (1V.27|) . we reproduce the abelian current conservation law (1V.15j) - (lV.21j) . 



The hydrodynamics of the AdS black brane is conformal, due to the conformal structure 
of the AdS boundary. We now wish to briefly discuss the subject of conformal symmetry, 
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which was not assumed in our derivation. For ideal hydrodynamics to be conformal, all that 
is required is a scale-covariant equation of state, i.e.: 

e(s,n a ) = s 4/3 e(n a /s). (VI.l) 

For the (leading-order) viscous hydrodynamics to be conformal, there are two requirements. 
First, the shear viscosity 77, the conductivity a a b and the coefficient £ a of the vorticity term 
in the charge current must likewise be scale-covariant functions: 

r) = sfj(n a /s); a ab = s^a^rf/s)- & = s 2 l%{n h /s). (VI.2) 

The second requirement is for the bulk viscosity to vanish: 

C = . (VI.3) 

In principle, we can have a system that is conformal in the ideal approximation, but loses 
conformal symmetry when the viscous corrections are taken into account. However, our 
results in the abelian case show that for the class of fluids described by an Einstein-Maxwell 
null horizon, conformal symmetry of the equation of state is sufficient for the viscous dynam- 
ics to be conformal as well. Indeed, eq. ( 1VI.3I) is satisfied always, while the explicit formulas 
for 77, a a t) and £ a guarantee that ( 1VI.1I) implies ( 1VI.2I) . In the non-abelian case, conformal 
symmetry at the viscous order does not quite follow from a conformal equation of state, 
since a^ b can be arbitrary. This will be the case, however, if the corrections to A a ^ which 
determine a^ b are governed by a conformally invariant condition, such as the condition ( II. 5j) 
used in AdS/CFT. 

In a previous version of this work, conformal invariance was assumed and used from the 



start. One can then use a Weyl-covariant formalism, as described in [28]. The main simpli- 
fication in that case is that the cancellations <jV758|> - < 1VT59|> (but not < N\57\i ) are no longer 
relevant: the Weyl-covariant equivalents of the corresponding terms vanish automatically, 
regardless of their coefficients. 

Finally, we would like to address the effect of gauge choice on our results. The horizon 
fields A® and j" are subject to 4d gauge transformations, with A ^ transforming as a connec- 
tion and j° transforming homogeneously. We are interested in descriptions where A^ and j° 
are local functionals of the hydrodynamic fields (h IJiU ,£ tJ- ,s,n a ). Therefore, we are only con- 
cerned with gauge transformation that preserve this property. These are the transformations 
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e lA for which the angle parameter A a is itself a local functional of (h fiu ,£ fl ,s,n a ). In gen- 
eral, the functional dependence of and j° on (h^, i^ 1 , s, n a ) will be altered by such gauge 
transformations. The hydrodynamic equations, however, remain unchanged at the first two 
orders we've been considering. To see this, note that the highest-order corrections relevant 
for these equations are A^ a and j^ a . Thus, we only need to consider gauge rotations with 
angles A a of order 1, e and e 2 . 

First, let us take A a ~ e 2 . To second-order accuracy, the only effect of such a transfor- 
mation is to rotate the zeroth-order fields A^ a and jl ^. Now, the second-order corrections 
appear in our derivation only in the combination f a bc{A^f >h j^ c ^ + A^' b j^ CIM ). The contri- 
bution to this quantity under the second-order gauge transformation should come from the 
rotation of f a b c A^ b j^ ctl = y/—hf a i, c fj, b n c , which vanishes. We conclude that second-order 
gauge transformations have no effect on the hydrodynamic equations. 

Let us now consider zeroth-order and first-order transformation parameters A a . To main- 
tain the first-order variable-fixing conditions of section IV B| the hydrodynamic fields must 
also transform under A a , in the trivial (h^, s) and adjoint (n a ) representations. Now, 
recall that all zeroth-order and first-order functionals A a necessarily lie in Z(n) (see sec- 
tion IV A II) . This means that the corresponding gauge transformations leave the fields 
(h flI/ ,£ l \s,n a ) unaltered. The hydrodynamic equations must therefore remain unchanged 
as wel]j. 
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3 With a little work, this can also be seen explicitly from the derivation of the equations. 
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Appendix A: Notations 

The black brane's event horizon (or a more general null hypersurface) is denoted by %. 
We use 5d coordinates x A = (r, where r = and d^r = 9ab^ B on H. The Yang-Mills 
gauge group is Q. 

1. Index conventions 

• Indices in the 5d spacetime's tensor bundle are denoted by uppercase Latin letters 
(A, B, . . .). They can be implicitly raised and lowered with the spacetime metric qab 
and its inverse g AB . 

• Indices in the 4d horizon's intrinsic tensor bundle are denoted by lowercase Greek 
letters (/x, u, ...). They are implicitly raised and lowered with the "hydrodynamic" 
metric h^ v and its inverse h^ u . Indices are never implicitly raised or lowered with the 
degenerate horizon metric 7^. 

• Indices in the adjoint representation of the gauge group are denoted by lowercase Latin 
letters from the beginning of the alphabet (a, b, . . . ). For non-Abelian components of 
the gauge group, an orthonormal basis with respect to the Killing metric is implied. 

• Lowercase Latin indices from the middle of the alphabet (i, j, . . . ) run from 1 to rank Q, 
and enumerate the independent scalars iVj and vectors (A^) a that can be constructed 
out of a gauge vector n a . 

2. Bulk quantities 

• e ABCDE is the 5d Levi-Civita density with components ±1. e abode is the correspond- 
ing inverse density, also with components ±1. 

• Qab is the 5d metric, with inverse g AB and determinant g. Rab is the corresponding 
Ricci tensor. 

• T^ u is the bulk stress-energy tensor. 
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• £ A is a vector field along the horizon's null generators. qab£ b is a covector tangent to 
the horizon. 

3. Horizon geometric quantities 

• e^ upa and e^ upa are the direct and inverse 4d Levi-Civita densities. 

• £ M is a vector field along the null generators. 

• S 11 — vi p is the area density current, independent of the scaling of l p . v is a scalar 
density that scales inversely with l^. 

• 7 MI/ is the horizon's degenerate metric. It is the pullback of qab into %. 

• = V M f is the Weingarten map. Its trace is n + 9. k is the surface gravity, and 9 is 
the expansion coefficient. The lowered-index version M P 7 P ^ = 9^ v = + (9/3)^^ 
is the horizon shear/expansion tensor, affi is the horizon shear tensor. 

• Qn" — ^(^V — K^a) i s the density whose divergence is used in the Gauss-Codazzi 
equation. Q^ v is its raised-index version with respect to 7^. 

• = \^^v — b p b Vl with arbitrary scalar and covector fields (A, b^), is an auxiliary non- 
degenerate metric which can replace 7 Ml , in certain formulas. Its inverse is (G~ l ) fiu . 

4. Yang-Mills quantities 

• A a A is the gauge potential. A a ^ is its pullback to "H. 

• fabc = f[abc] are the gauge group's structure constants, with the coupling strengths 
included. 

• Pabc — P(abc) are the Chern-Simons coefficients. They form a gauge-invariant constant 
tensor. 

• F a AB = 2d [A A a B] + f a hc A h A A% is the gauge field. F° v = 2d {il A a u] + S a ^K is its 
pullback to the horizon. jg = {yJ—g/A^F^ is a 4d vector density on "H, related to 
the hydrodynamic charge current. 
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5. Charged brane/hydrodynamic quantities 



• m and q a are mass and charge parameters of the black brane. 

• is the 4d metric of the hydrodynamics, with inverse h^ u and determinant h. It can 
be either flat or curved. 

• = yZ—hsi^ is the entropy current, is the entropy velocity, normalized with respect 
to h^ u . s is the entropy density in its rest frame. £ M is the same as the horizon's null 
generator defined in section IA 31 Note that = h^f is not the pullback of gAB^ B to 
% (which is zero). By the entropy-area relation, we have s M = S^/A. 

• = 5% + l^lv is the projector orthogonal to ^ (with respect to h^J). 

• R is the brane's horizon radius, in the sense that s = R 3 /4. 

• J a/i is the conserved hydrodynamic current. n a ~ q a is the charge density in its rest 
frame. J^g is the axial part of J aM , derived from the Chern-Simons term. a a b is the 
conductivity matrix in the Landau frame. a a b is the conductivity matrix in the entropy 
frame. 

• is the hydrodynamic stress-energy density, conserved with respect to h^. is the 
fluid's energy velocity (timelike eigenvector of Tjf), normalized with respect to h^. 
e ~ m is the mass/energy density in its rest frame, p is the pressure. 

• T is the temperature of the brane/fluid. fi a is the chemical potential. 

• Z(n) is the centralizer of n a , i.e. the subspace of the gauge algebra that commutes with 
it. Z a i, is the projector into Z(n). Z^ b is the projector into its orthogonal complement 
Z». 

• Ni is an exhaustive set of charge scalars constructed out of n a . The charge vectors 
(iV;) a = dNi/dn a constitute a basis for Z(n). 

6. Derivatives 

• (or Oa) is the ordinary partial derivative. Since we'll be dealing with distinct 
covariant derivatives, we use the partial-derivative symbol to highlight connection- 
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independent objects. 

Lt is the Lie derivative with respect to P. 

V a is the bulk covariant derivative, associated with the gravitational metric #ab- 
is the restriction of into "H, in contexts where it is well-defined. 

D p is the covariant derivative associated with the 4d "hydrodynamic" metric h^. The 
derivative along £ M is denoted by D = £^D p . 

is the fluid's expansion rate. Di p is the acceleration, = P?P° D^ p £ a ) is the 
fluid's shear tensor, uo pv = PPP^D^l^ is the vorticity. = (l/2)e^ vpa i u uo p(T is a 
vector density which also describes the vorticity. To the relevant order, these are all 
equivalent to the analogous quantities defined in terms of the energy velocity vP 1 . 

Djj, is the covariant derivative with respect to the auxiliary metric G pu . 
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